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We study a classical double exchange magnet with direct antiferromagnetic superexchange coupling, J, 
between the localized spins. It is shown that the de-stabilization of the ferromagnetic ground state with 
increasing J leads to phase separation; the latter always preempts the spin-wave instability (softening of the 
magnon spectrum) . It is also found that the boundaries separating the ferromagnetic and antiferromagnetic 
areas of the sample tend to be abrupt. 
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The phenomenon of phase separation in double ex- 
change (or s-d exchange) magnets has been discussed in 
the literature for a long time [1]. Presently it attracts 
attention of both theorists and experimentalists investi- 
gating the properties of doped manganese oxides which 
exhibit colossal magnetoresistance (CMR) [2]. 

One of the simplest yet physically relevant scenarios 
of phase separation takes place in a single-orbital double 
exchange magnet in the presence of a sufficiently strong 
antiferromagnetic superexchange coupling J between the 
neighbouring localized spins. Some of the earlier arti- 
cles [3-5] noted that the thermodynamic instability of de 
Gennes' canted phase (originally [6] thought to emerge 
as an outcome of double exchange - superexchange com- 
petition) at small values of band filling x indicates that 
this phase is unstable with respect to phase separation. 
Others [7,8] used mean field and variational approaches 
to minimise the ground state energy and obtained rich 
phase diagrams. The numerical studies [2,9] suggest that 
phase separation occurs at all values of a;, provided that 
the value of J is sufficiently large. The systematic physi- 
cal understanding of phase separation in these systems is 
however still incomplete. In the present article, we will 
address some aspects of this deficiency, paying special at- 
tention to the relationship between phase separation and 
spin waves, as well as to the features of the interphase 
boundaries, which are likely to affect both equilibrium 
and transport properties of the system. 

We consider a model defined by the Hamiltonian 
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Here Cja (with a =t,i) are the electron annihilation 
operators, Si are the operators of the core (localized) 
spins located at the sites of a square lattice, and the 
vector <7°-^ is composed of Pauli matrices. We assume 
that ionic spins are classical [S 00), and the Hund's 
rule exchange coupling is large {Jh — 00). Under these 
assumptions, there is only one carrier spin state avail- 
able at each site, and the hopping coefhcient is effectively 
renormalized [6,10] by a factor [S^ + 5, • Sj)'^/^ /{SV^). 
Throughout the paper we use units in which hopping t 
and the lattice spacing are equal to unity, and we con- 
sider the r = case. While technically our treatment is 
restricted to the two-dimensional model, we expect that 
the conclusions remain qualitatively valid in three dimen- 
sions. 

When the value of superexchange J is small, the 
ground state of the system is ferromagnetic (FM), and 
the carrier spectrum is given by [11] ej: — — cosfc^; ~ 
cosfcy. The total energy of the conduction band, E = 
(where ^(e) is the corresponding density of 
states) is determined by the value of chemical potential 

As the value of J increases, the ferromagnetic ground 
state of a double exchange magnet is eventually destabi- 
lized, which can be seen in particular from the fact that 
at a certain J^^ = \E\/8 the spin wave spectrum of FM 
phase LUj} = (-2J + \E\/4){2 + Cjij/S (see Ref. [12] and 
references therein), softens and then turns negative. The 
nature of emerging phase has not been clarified by a spin 
wave theory investigation [12], although it was demon- 
strated that the homogeneous de Gennes' canted state 
can never be stabilized. We will now show that in re- 
ality this spin wave instability of a homogeneous, fully 
saturated FM state can never be reached as the system 
becomes unstable with respect to phase separation al- 
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ready at J < Jsw This conclusion will be verified for 
all values of band filling x, and the emerging new phases 
will be discussed. 

As the value of J increases, the homogeneous FM 
state with a chemical potential ji = i^FAiix) is desta- 
bilized with respect to phase separation whenever its 
thermodynamic potential, QpM = E — /ix + 2 J, be- 
comes larger than the thermodynamic potential of some 
other phase, calculated at the same value of /U. Fur- 
ther increase of J will result in an increase of the vol- 
ume occupied by this new phase at the expense of the 
FM area; the combined change of values of J and ^ will 
eventually lead to another instability or phase transi- 
tion, etc. Thus in order to determine the stability re- 
gion for the homogeneous FM state at a given value of 
X, one has to identify the phase P, for which the condi- 
tion flpd^pM^x), Jp) = ^PM{fJ'FM{x), Jp) is reached at 
the lowest possible Jp. In doing so, we are restricted by 
the usual limitations of variational approach [13]. Re- 
sults are summarized in Fig. 1; note that for any [14] 
value of X, there is at least one phase corresponding to a 
value of Jp smaller than Jsw{x). This allows us to con- 
clude that phase separation, rather than a second-order 
phase transition associated with the spin- wave instability 
at J = Jsw , is the generic outcome of double exchange - 
superexchange competition. 

We consider the following phases: 
G-type AFM (G-AFM) phase with checkerboard spin 
ordering. The thermodynamic potential is given by [15] 
Q,G = — 2J. Given our selection of phases, the G-AFM 
phase corresponds to the lowest value of Jp at < a; < 
0.245. 

A-type AFM (A-AFM) phase characterized by the 
wave vector {7r,0}. We find that this phase with ther- 
modynamic potential fi^ = — {^/l — H'^ + yitarcsin/i + 
7r/i/2}^(l — /x^)/7r yields the lowest value [16] of Jp at 
0.37 <x < 0.5. 

Island phases represent a generalization of one- 
dimensional island phases [17], and were already found 
in the numerical investigations of both two- and three- 
dimensional cases [9]. The thermodynamic potential 
for the 2x2 island phase (see Fig. 1 6) is given by 
^2 = —{l^ + 1)^(1 — /^^)/4, and this phase corresponds 
to the lowest Jp{x) for 0.293 < a; < 0.335. Similarly, at 
0.335 < x < 0.37 the lowest value of Jp corresponds to 
the 3x3 island phase. 

Chain phase (see Fig. 1 h) is proposed here for the first 
time. Like in A-AFM phase, carriers in the chain phase 
arc not completely localized: they can travel along the 
braided paths parallel to the direction of the chains. We 
find 

f^c/i = |\/2-/x2 + ^arcsin^-f |/^| 6{2~^i'')+^- , 

and the corresponding value of Jp is the lowest for 
0.245 <x< 0.293. 



Phase separation obviously leads to formation of inter- 
phase boundaries separating different areas of the sample. 
It is very unlikely that these are smooth, Bloch-like walls 
with gradually increasing magnetization and decreasing 
staggered magnetization - indeed, this would inevitably 
involve the appearance of (unstable [12]) canted spin or- 
dering in the intermediate region. The extreme opposite 
case of an abrupt interphase boundary corresponds to a 
complete change of magnetic order over one lattice link 
(see Fig. 2 a); for the case of phase separation into the A- 
AFM and FM phases, such boundaries should run along 
a crystal axis, and for the G-AFM or chain phases - in a 
diagonal direction. Short of a rigorous proof that these 
abrupt walls represent the optimal configuration, we will 
probe their stability with respect to small "smearing" 
perturbations [18] shown schematically in Fig. 2 a. The 
energy cost (per unit length) of an abrupt wall separat- 
ing the FM phase form another phase P can be written 
in the form 

5f2 = Wp{n, J) + Ap{iJL, J)^^ + Bpi)'^ , tp, ^ « 1 . (2) 

The evaluation of W, A, and B is sketched in the Ap- 
pendix. 

The values of the coefRcients Ap{iJiFM(x), Jp{x)) and 
Bp{ijlfm{x), Jp{x)) are plotted in Fig. 2 h. We see that 
the energy cost of perturbing an abrupt FM G-AFM 
boundary is always positive throughout the region where 
the value Jg{x) is the lowest among the Jp's evaluated 
above (sec Fig. 1 a). A similar situation occurs for an FM 
- A-AFM boundary. This represents rather convincing 
evidence to the effect that the boundaries separating the 
FM and AFM areas are abrupt. This conclusion is not 
altogether unexpected as it is known, e.g., that small 
polarons propagating in a doped antiferromagnet have 
no tails [19], i.e., arc "abrupt". 

We speculate that the instability of an abrupt chain - 
FM wall at a; > 0.265, seen from Fig. 2 b, may indicate 
the existence of another, so far unidentified phase P with 
a lower critical value of superexchange, Jp < J chain- 
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APPENDIX A 

Here we outline the main steps in the calculation of 
coefficients on the r. h. s. of Eq.(2). Since it is assumed 
that the interphase boundary is perfectly straight, the 
variables separate, and one is left with an effective one- 
dimensional problem (Fig. 2 a). For simplicity, we will 
describe the case of a boundary between the FM and G- 
AFM or A-AFM phase; the chain phase is treated simi- 
larly. Carrier hopping in the direction perpendicular to 
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the boundary vanishes in the AFM phases, so the pres- 
ence of a bulk AFM phase to the left of the boundary does 
not affect the electron wave functions in the FM part of 
the sample. The electronic contribution to the energy 
cost of a boundary is due solely to the associated change 
in the wave functions in the FM part of the sample (these 
do not extend across the boimdary). Thus, the energy of 
the interphase boundary is closely related to the energy 
of an abrupt FM domain wall depicted in Fig. 2 c (top); 
indeed, it is obvious that the latter energy is equal to 
2Wp + Apifi^ + Bp^^. The origin of the prefactor 2 in 
the first term is clarified in Fig. 2 c (bottom), showing 
two FM-AFM boundaries (only one of which is perturbed 
at (p,ip 0), separated by three simple chains of AFM 
between them. Generally, the energy difference between 
such "stripe" configurations and the domain wall shown 
in Fig. 2 c (top) is equal to ^Iafm — ^fm times the 
number of inserted chains of AFM [20], i.e. is entirely 
due to the insertion of the bulk AFM phase. In turn, 
the abrupt domain wall is a local perturbation of the FM 
order, and its energy can be conveniently evaluated with 
the help of the Lifshits-Krein trace formalism (see Ref. 
[4], and references therein). We find [21]: 
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and rather lengthy expressions for Ap and Bp. At 3; <C 1, 
we obtain Wq ~ -V2J + {n + 2)3/2/(3v^7r); the 3D 
counterpart of this expression was used in Ref. [1] . 
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Misprints in the published version of these equations are 
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FIG. 1 . (a) Values of superexchango J corresponding to the 
instabilities of FM order. The soUd Une labeled "G" ("A") 
represents the phase separation into G-type (A-type) AFM 
and FM phases, whereas the dashed line labeled "2" ("3") 
corresponds to phase separation into 2x2 (3x3) island phase 
and FM phase. Phase separation into FM and chain phases 
(dashed-dotted line) may also be possible at certain values of 
band filling x. Note, however, that the spin-wave instability 
at J = Jsw(x) (dotted line) is always circumvented by phase 
separation, (h) Spin ordering in the 2x2 island (left) and chain 
(right) phases. 
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FIG. 2. (a). Spin arrangements near the abrupt bound- 
aries between FM and G-AFM or A-AFM phases (top) and 
between FM and chain phases (bottom). Each arrow repre- 
sents a FM chain parallel to the boundary; interchain dis- 
tance is equal to unity for the FM - A-AFM boundary, and 
to l/v/2 for the FM - G-AFM and FM - chain cases. Dashed 
arrows correspond to perturbed boundaries with <p, 1/) 7^ 0. 
(h) The values of coefficients Ap{iifm{x), Jp{x)) (left) and 
Bp{hfm{x), Jp{x)) (right) for the boundaries between FM 
and G-AFM (solid), A-AFM (dashed), or chain (dotted) 
phases. See Eq. (2). (c) An abrupt domain wall in the FM 
state - see Appendix. 
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